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Abstract 

> . 

The enhanced suppression of J / ifi production at large xp in pA collisions is studied 
in the framework of gluon depletion at large x\. The nonperturbative process that 
^f) ' modifies the gluon distribution as the gluons propagate in nuclear matter is described 

by an evolution equation with a kernal to be determined by phenomenology. With 
nuclear shadowing and anti-shadowing taken into account, the effect on the gluon 
distribution is shown to be a depletion in excess of 40% at x\ ~ 0.8 for A > 100. There 
is a small amount of enhancement of the gluon distribution at small x\, but it does 
not lead to any contradiction with the existing data on J/ip suppression in the central 
region. Extentions to V' suppression and AB collisions are also investigated in the 
framework of gluon redistribution. 



^ : 1 Introduction 

In an earlier paper |l] we presented the phenomenological evidence for the depletion of high- 
momentum gluons as a projectile proton traverses a target nucleus in pA collisions. Here 
we present a more complete discussion of the subject with more technical details and with 
extension to nucleus- nucleus collisions. 

Our approach is unconventional in that we do not make the usual assumption that the 
parton distributions in a proton remain unaltered as the proton propagates through a nu- 
cleus, even if a hard subprocess occurs deep in the nucleus. That assumption amounts to 
factorization, a property that has been proven for pp collision, but clearly must fail for pA 
collision if A is infinitely large. For realistic nuclear sizes our alternative assumption that the 
parton distribution can be modified led to J/ip suppression that cannot be distinguished 
from the effects of the usual mechanisms |3], What is different now is the appearance of 
new data on J/ip suppression in pA collisions at large xp 0. Since those data cannot be 
explained in terms of hadronic absorption of the produced cc state || 0, §], we find more 
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direct support for the idea of gluon depletion before the hard subprocess that produces the 
cc state [|TJ. 

The data of the Fermilab E866 experiment || on the J/ip suppression in pA collisions at 
800 GeVc are expressed as the ratio 

R (x F , A) = a A (x F ) /Aa N (x F ) = A*^ 1 , (1) 

where <Jn,a is the cross section for J/ip production by a proton on a nucleon (N) or on a 
nucleus (A). An analytical formula for a(x F ) is given in || 

a(x F ) = 0.96 (l -0.0519^- 0.338:4) , (2) 

for —0.1 < x F < 0.9. Equation (^) differs significantly from the one given in their preprint 
||, on the basis of which the analysis in [|IJ was done. Here we perform a reanalysis based 
on the new parameterization in Eq. (Q). 

In our view there are three sources that can contribute to the x F dependence, and we 
write them in a product form 

R(x F ,A) = G(x F ,A)N(x F ,A)H (x F ,A) , (3) 

where H(x F ,A) represents the hadronic absorption of the cc state, N(x F ,A) nuclear shad- 
owing, and G(x F , A) the gluon depletion effect, defined as 

G(x F (xi),A) = g(x h A)/g(xi,0) , (4) 

X\ being the momentum fraction of the gluon in the projectile proton, and g{x\, A) being the 
effective gluon distributions of that projectile in a nucleus A at the point of cc production, 
averaged over the penetration depth at which the hard process occurs. 

By H(x F , A) we mean the absorption factor that operates between the production of the 
cc state and the detection of the final J / ip due to any mechanism, including the interaction 
with comovers. Since no good arguments have been advanced to show that H(x F ,A) can 
have a significant dependence on x F || 0, U, we shall assume in the following that H(x F , A) 
is independent of x F . This is not a serious limitation in our formalism. If a reliable x F 
dependence is found at a later date, we can easily incorporate it in our analysis. For now we 
adopt the usual Gerschal-Hufner form. 

H(A) = exp [-apL(A)} (5) 

where a is the absorption cross section, p the nuclear density, and L(A) the mean path length 
in A that a cc state propagates. 

In the next section the nuclear shadowing factor N(x F , A) will be discussed. We shall find 
a simple formula that can represent the change in the gluon distribution in the target nucleus 
due to shadowing and anti-shadowing. Such a formula facilitates the analysis, and offers a 
simple parametrization that is convenient to use, independent of the particular problem that 
we apply it to here. 

The determination of the depletion factor G(x F (xi), A) is the main theme of this paper. 
In Sec. 3 we shall go beyond a review of the content of Ref. |jj , not only because the data has 
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changed from those in the preprint 0, resulting in numerical differences, but also because 
we shall improve on |]] in some technical details and include some new material. 

In Sec. 4 we shall extend the result from the study of the pA problem to AB collisions. 
We shall show how the enhanced depletion at large xf does not affect the J/ip suppression 
at small Xp, which is where the existing data for AB collisions were collected. We shall have 
nothing to add on the subject of enhanced suppression observed in Pb-Pb collisions beyond 
what we have advanced in Ref. @. 



2 Nuclear Shadowing and Anti-shadowing 

The nuclear shadowing and anti-shadowing problem has been studied phenomenologically 



by Eskola et. al. [10, |TT|. Instead of focusing on the physics of the origin of the problem in 
QCD, they analyzed the deep inelastic scattering data of nuclear targets. On the basis of 
DGLAP evolution |12| they can determine the parton distribution at any Q 2 > 2.25 GeV 2 . 



The results are given in terms of numerical parametrizations (called EKS98 [fLl]| ) of the ratio 

Nf(x,Q 2 ) = fa/ A (x,Q 2 )/fa(x,Q 2 ), (6) 

where fa is the parton distribution of flavor i for the free proton and fa/ a is that for a proton 
in a nucleus A. 

For the purpose of J/i/j production we are interested in Eq. for only the gluons. We 
consider only the dominant subprocess g\{xi) + g 2 (x 2 ) ~ ¥ c + c, where X\ is the momentum 
fraction of the projectile gluon whose distribution is gi(xi), and x 2 is that of the target gluon 
whose distribution is g 2 {x 2 )- For the cc state produced with momentum fraction xp = X\ —x 2 , 
the usual kinematical relations are 

x h2 = [(Ar + x F ) 1 / 2 ±x F }/2, x 1 x 2 = t = M 2 j/iP /s, (7) 

where it is assumed that the cc state that turns to J/ip is produced near threshold. Thus 
the virtuality of the subprocess g + g c + c is given by the value of Mj/^, or slightly 
higher. We shall take Q 2 = 10 GeV 2 , a value that is chosen in EKS98 ]TT] to give explicit 



values of N i A {x\ 1 Q 2 ). For simplicity we shall label the gluon distribution at Q 2 = 10 GeV 2 
by N(x 2 ,A). 

The numerical output of EKS98 for N(x 2 , A) is shown by the points in Fig. 1 for A = 50, 
100, and 200. We exhibit only the values for x 2 in the range 0.01 < x 2 < 0.12, since that is 
the range relevant for the production of J/ip at 800 GeV/c for < x F < 0.8. For the purpose 
of convenience in using those values of N(x 2 , A) in analytic manipulation and computation, 
we propose a simple formula that contains the shadowing and anti-shadowing effects. Since 
the cross-over of the two effects occurs at x 2 = 0.02 where N(x 2 ,A) = 1 for all A, it is 
sensible to use an auxiliary variable £, defined by 

f = 3.912 + lnx 2 (8) 

which vanishes at x 2 = 0.02. Moreover, we notice that lniV(a;2,^4) depends linearly on InA 
to a good approximation, so it suggests a power-law behavior 

N{x 2 ,A) = A^\ (9) 
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The exponent (3(x2) can be determined by fitting the data for A = 100. In terms of £ we 
find a good fit with the parametrization 

/?(£) = £(0.0284 + 0.0008£ - 0.0041£ 2 ), (10) 

the result of which is shown in Fig. 2, where the points are for hiN(x2,A)/hiA obtained 
from EKS98 at A = 100. This convenient formula for /?(£) can then be used in conjunction 
with Eq. (H) to determine N(x2,A) for other values of A. The curves in Fig. 1 exhibit the 
good agreement between the data and our parameterization for A = 50, 100, and 200. 

In the following we shall simply use Eq. (|j) as a summary of the effects of nuclear 
shadowing and ant i- shadowing for problems in pA collisions where the relevant range of x<i 
is in the interval 0.01 < Xi < 0.12. 



3 Evolution of Gluon Distribution in a Nucleus 

We now may regard R(xp,A) and N(x2,A) as known phenomenologically. Thus, from Eq. 
(|) we may write 

G(x F , A)H(A) = A«(^)-/3(*2(^))-i. (H) 

Although the form of H(A) is given by Eq. (|5|), we do not know the value of o in the 
present circumstance where we allow the possibility of gluon depletion. Hence, we treat 
H(A) temporarily as unknown, along with the depletion factor G(xp,A). However, the Xp 
dependence is completely known from the RHS of Eq. (|TTD . 

To proceed we need some theoretical input on the possible form of G(xp,A), or more 
directly the gluon distribution g(x±, z) in the projectile, where z is the distance traversed 
in a nucleus. Note that we have refrained from referring to the projectile as the proton, 
since the possible modification of g(x\,0) for z > implies that the incident proton loses 
its usual identity, in particular, the nature of its partonic content, as the projectile, now 
identified only as a flux of partons, propagates in the nuclear medium. How g(xi, z) evolves 
in the nuclear medium is clearly a nonperturbative process that involves multiple scatterings 
of gluons and quarks at low virtualities. Nevertheless, for every incremental distance, dz, 
that a gluon travels the modification that g(xi,z) undergoes must be perturbative in that 
g(xi, z+dz)—g(xi, z) is small and is proportional to dz. It is therefore reasonable to adopt an 
evolution equation similar in spirit to that of DGLAP [T2|], but with the change in resolution 



scale dhiQ 2 replaced by the change in penetration depth dz, so that we write 

4-9(*,*)= f 1 ^g(x',z)Q(^), (12) 
dz Jo x' x' 

where the unknown kernel Q(x/x') controls the gain and loss of the gluons in dz. Q{y) 
cannot be determined by perturbative calculation, as the splitting functions in pQCD for 
Q 2 evolution. Equation ( |T2] ) is similar to the nucleonic evolution equation proposed in |I3 



except that this is now at the parton level. In Eq. (|T2| ) the quark sector has been left out for 
simplicity. To be more complete one should include also the effects of the couplings of gluons 
with the quarks, a task that is deferred to the future. Thus what we can achieve now is the 
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determination of an effective kernel Q(y) that can account for the enhanced suppression of 
J /if) at large xp. 

To determine the z dependence of g(x, z) let us take the moments by defining 

9n{z) = [ dxx n ~ 2 g(x,z) (13) 

and 
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Q n = f 1 dxx n - 2 Q(y). (14) 
Jo 

Then by the convolution theorem, Eq. (|T2|) becomes 

dg n (z)/dz = g n (z)Q n , (15) 

whose solution is 

g n (z) =g n (0)e zQ ". (16) 

It is possible that Eq. fll5|) is valid only when z is large enough, in which case Eq. (|l~6|) should 
be modified to read 

9n(z) = g n (z )e^ Q - (17) 

for z greater than some positive value of Zq. 

To proceed, let us substitute Eq. (f|) in flTl"| ) and define 

J(xi, A) = g(xi, o)A a(xF(xi)) - /3(x2(xi)) - 1 (18) 

where the interrelationships among Xi,x 2 and xp are specified by Eq. (0). For g(xi,0) we 
use the canonical form 

g(x u 0)=g (l- Xl ) 5 . (19) 

The final result is insensitive to its form and independent of go, which we shall set to be 1. 
Thus we may regard J{x\, A) as known. Since we also have 

J(x 1 ,A)=g(x 1 ,A)H(A), (20) 

its moments are, by virtue of Eqs. (f|) and (|16|), 

J n (z) = g n (0)exp[z(Q n - op)]. (21) 

Here and in the following we shall use z (until Sec. 5) to denote the average penetration 
depth (i.e., z = za) in A when a cc state is produced by gg annihilation. It is then also 
the average length that the cc must travel in A and be subject to hadronic absorption, i.e., 
z = L(A). 

We can determine J n (z) by taking the moments of J(xi,A), as expressed in Eq. (|18D . 
However, there is a problem in evaluating J n (z) = dxiXi n ~ 2 J{x\, z), since J(x\,A) is 
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Table 1: Values of the coefficients ki 



A 


z 


k 


ki 


k 2 


h 


A4 


100 


4.177 


-0.5444 


-0.2063 


5.87 x 10- 3 


-8.35 x 10- 5 


0.424 


200 


5.262 


-0.6267 


-0.2354 


6.64 x 10~ 3 


-9.39 x 10- 5 


0.486 



ill-defined at x\ = 0. According to Eq. (|7|), x 2 diverges as x\ — > 0; thus limiting x 2 to 1 
implies that x\ cannot be less than Mj,^/s, a small but nonvanishing value. Furthermore, 
xp becomes negative at small x\, and we lose any knowledge about a(xp) for xp < —0.1 [|J. 
Also, (3{x 2 ) is not reliably known at large x 2 , so the RHS of Eq. ( |T8| ) cannot offer accurate 
determination of J(xi,A) as x% — > 0. These defects can be suppressed by the factor Xi n ~ 2 
in the integrand, if we restrict n to > 3. We shall therefore determine J n only for n > 3. 
For convenience, let us define 

K n {z) = ln[J n {z)/g n (0)} = z{Q n - op). (22) 



Using Eqs. (0), (|T0D and (jl9f) in ([18|), we can calculate J n (x) and therefore K n [z) for n > 3. 
As mentioned earlier, Eq. is different from a previous form of a(xp) given in M and 
used in 0. The results are shown as discrete points in Fig. 3 for A = 100 and 200. The 
corresponding values of z are halves of the average total path lengths of the nuclei, i.e., 
z = 3Ra/4: = 0.9A 1 / 3 fm and are therefore z\ = 4.177 and z 2 = 5.262 fm, respectively. 

To extract the information contained in those points in Fig. 3 we need an analytical 
representation of K n . We choose to fit K n by the following formula, different from the one 
used in [0, 

3 

K n = J2k i n i + k 4 n 1/2 . (23) 

i=0 

The results of our fits are shown by the smooth curves in Fig. 3. The corresponding param- 
eters are given in Table 1. Note that the fits allow us to extrapolate smoothly to n — 2, 
where we could not calculate J 2 . 



Before we determine Q n from K n in Eq. (|22), we need to verify the z dependence. It 



should first be recognized that the experimental parameterization of the A dependence, such 
as in Eq. (JI|), is not compatible with the theoretical expectation, such as in Eq. (||) and (|2T|), 
except in a certain range of A. Since a power-law A 7 , expressed as e 7lnA , is approximately 
e 7 ' 2 , where z oc A 1 / 3 , only when \nA is approximately A 1//3 , the correspondence can only 
be for 60 < A < 240. With that understanding, let us nevertheless calculate J(xi,z) for 
all z < 6 using A = (z/0.9) 3 in Eq. flT8p, take the moments, and then determine K n (z) 
through the first half of Eq. (p2|). The result is shown as points in Fig. 4 for eleven values of 
z between 0.9 and 5.9, corresponding to A from 1 to 282, and for four representative values 
of n, viz., 3, 8, 13 and 20. The straight lines are linear fits of the last six points for each 
value of n. Evidently, the z dependence of K n (z) is very nearly linear for 3.4 < z < 5.6, 
which corresponds to 54 < A < 240. Thus our theoretical formalism is consistent with the 
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experimental data in the region where In A ~ A 1 ^ 3 . At z = 0.9, or A = 1, all points converge 
to K n = 0, as they should. We cannot reliably apply our formalism to the collision problems 
where A < 50. Fig. 4 also suggests that even when A is large, say > 100, the gluon evolution 
equation (|l^) may not be valid at small z, here used in the sense of penetration depth within 
the large nucleus, not the average depth. In the following we shall limit our consideration to 
only the linear portion of Fig. 4. In that region the second half of Eq. ( p2|) should be treated 
as differentially correct, i.e., 

AK n {z)/Az = Q n - pa (24) 

where AK n (z) = K n (z + Az) - K n (z). 

Since the values of ki in Table 1 are determined in the linear region, we can use them to 
obtain Q n . If we write 

3 

Qn = E^ + ^ 1/2 , (25) 
i=0 

then we have from Eqs. (p3|) and ([24]) 

g = Ak /Az + pa, q t = AkJ Az (1 < i < 4). (26) 

For the two z values in Table 1, we get (with Az = 1.086 fm) 

q = -0.0758 + pa, q x = -0.0268, , . 

q 2 = 7.13 x 10- 4 , q 3 = -9.6 x 10~ 6 , g 4 = 0.0568 1 ' 

in units of fm" 1 . 

Since the absorption cross section a is unknown when gluon depletion is not negligible, 
q Q is not fixed by Eq. (p7[). Whatever the dynamics of gluon depletion is, we require that the 
total gluon momentum does not increase with z. Since the gluon momentum is / dxg(x, z) = 
g2(z), that requirement implies in conjunction with Eq. ([17]) that Q2 < 0. Choosing the upper 
bound Q2 = leads to the condition, on account of Eqs. (|25|) and 



-(2gi + 4g 2 + 8g 3 + V2q A ) = -0.0295. (28) 



With all parameters in Eq. (|25|) now determined, the n dependence of Q n can be exhibited, 
as shown in Fig. 5. Evidently, Q n is smoothly varying and the use of the polynomials in Eqs. 
(p3|) and (|25"D, which are different from those used in ||], is justified. 

Although it is more direct to proceed immediately to the use of Eq. (0) to the deter- 
mination of g n {z) and therefore G(x\,z), which is our goal, there is some advantage in an 
attempt to find Q(y) at this point, while we are on the subject of Q n . The easiest way to 
do that is to put Q n in form 

Q n = c + J2c j /(n + j-l) , (29) 

j 

used in M, so that it can imply directly 

Q(y) = c 6{l-y)+ y £c j y i > ( 30 ) 
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To translate from Eq. ( p5|) and (p9"D, we fit the values of Q n at the 19 integer points, 2 < 
n < 20, determined by (|25|) , by use of the formula ( p9"D with a suitable number of terms in 
the sum. It turns out that a good fit can be achieved with three terms: j — 3,4 and 5. The 
result is 

c = -0.1988, c 3 = 6.205, c 4 = -23.316, c 5 = 19.866 . (31) 

We show in Fig. 6 both the discrete values of Q n at integral n and the fitted curve using 
Eqs. ( |29"D and (|31~1) . The corresponding function Q(y), calculated using Eqs. (^) and fl3~l|) , is 
shown in Fig. 7. It is evident that the Cq and C4 terms correspond to gluon depletion, while 
the C3 and C5 terms correspond to gluon regeneration. 

To determine G(x±, z) as defined in Eq. (|]), we have Eq. ( p!7[ ) that specifies the evolution 
in 2 in the linear region (see Fig. 4) from zq. At this point we have no formalism to extrapolate 
in the nonlinear region from z down to 0, and it is in reference to g(xi,0) that G(xi,z) 
is defined. However, in view of the fact that the hadronic absorption term H(A) is known 
empirically to be an exponential, as in Eq. @, i.e., \nH{z) being linear in z for all z, one may 
regard the nonlinear portion of Fig. 4 to be due primarily to the mismatch between \nA and 
A 1 / 3 at low A. Then we adopt the approximation that Eq. ([TB] ) is adequate for relating g n (z) 
to g„(0). With Q n being known from Eqs. fl2"5D, ( ]2"7| ) and (^), and g n {0) = B{n — 1, 6) which 



is the beta function, we can calculate g n (z). The result is shown in Fig. 8 by the full (open) 
circles for A = 100 (200), It is natural to fit the resultant g n (z) by a linear combination of 
beta functions in the form 

3 

g n (z)= y £a i {z)B{n-l,5 + i) , (32) 

i=l 

The coefficients ai(z) are determined by fitting g n {z)/g n {0) in order to reduce the range of 
variation. For the two values Z\ and z 2 , corresponding to A = 100 and 200, we obtain 

ai {z x ) = 0.3526, 02(^0 = 1.44, a 3 (20 = -0.78, , . 

01(23) = 0.2362, o 2 (2 2 ) = 1.655, a 3 (z 2 ) = -0.869. 1 j 

The curves in Fig. 8 are generated using Eqs. (|32|) and (^). Evidently, the fits are good. 
The inverse transform of the moments in Eq. (B2I) is 



g(x l ,z) = j2a l (z)(l-x 1 )^ , (34) 



i=i 



whose implication for 

G{x u z) = g(xi,z)/g(x 1 ,0) 



ai(z) + a 2 (-2)(l - xi) + a3(2)(l - xi) (35) 



can readily be calculated using Eq. Q3"3]). The results for z\ and z 2 are shown in Fig. 9. 
Clearly, there is significant depletion of gluons at large Xi, roughly 40% at X\ ~ 0.8. There is 
a small amount of regeneration at small X\. The cross-over is at around X\ ~ 0.2. Although 
the enhancement at small x±, is at the 2 to 3% level, in terms of the number of gluons in a 
small dx interval it is not insignificant compared to the depletion at large X\, because g(x±, 0) 
is strongly damped at large x\. That is how the condition Q 2 = is satisfied. 
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4 Suppression of J/ip and ip' 



Having determined G(xi,A) in the previous section, we can now return to the problem of 
charmonium suppression, including if)' . The Fermilab E866 experiment || gives data for 
both J/ip and if}' in the form of a{xp). In using Eq. (|3|) to calculate R(xf,A), and then 
a(xp), we have Eqs. fl35p , @ and (H) for G(x\,A), N(x2,A) and H(A), respectively. For 
both J/ip and ip' we use a = 6.5 mb in H(A). The results on a(xp) are shown as curves in 
Fig. 10 for J/ip and Fig. 11 for if}'. The two curves in each figure are for A = 100 and 200; 
they are sufficiently close to each other to be almost independent of A, thereby affirming our 
gluon evolution model for that range of A. 

The most significant part of what we have learned from this work on the pA data is that 
it is hard to reproduce the strong damping of the measured a(xp) at large xf without a 
substantial amount of gluon depletion at large x±, as seen in Fig. 9. Furthermore, if there is 
significant modification of the parton distribution as the penetration depth increases, then 
it is hard to justify the notion that a proton can traverse a large part of the target nucleus 
without any changes, or that it can be wounded upon the first collision with a nucleon and 
then remaining unchanged thereafter. 

So far out attention has been given to the xp > region only where data exist. However, 
the xp < region has interesting physics also, and the data can be obtained either at RHIC, 
or at a fixed target experiment with proton being the target for a heavy-ion beam. At 
large negative xp, x± would be small where Fig. 9 shows a very small enhancement, not 
depletion. On the other hand, X2 would be large, where N(x2,A) would exhibit not only 
anti-shadowing, but also the EMC and Fermi motion effects. As an illustration of the xp 
dependence, we have performed the calculation for J j ip suppression in pW collision at three 
energies and for all xp, negative as well as positive. The absorption cross section is set at 
a = 8 mb, a value found in the next section. The values of N(x2, A) for the full X2 range are 
obtained from EKS98 ||11|| . The results are shown in Fig. 12. We see that the negative xp 
region shows very little dependence on energy and reflects mainly the property of N(x2, A). 

Gluon depletion of the type discussed in this paper has no effect on the suppression of 
J /if) in the xp < region. However, there is another type of gluon depletion, discussed 
in the second paper in Ref. 0], called nonlinear depletion, that can influence the survival 
probability in the xp < region. When X2 is large, the fast gluons in the rear part of the 
nucleus can catch up and interact with the slow gluons released by the p-N collisions in 
the front part of the nucleus, leading to a depletion of gluon at large x 2 . Although this 
involves the interaction between gluons in different nucleons in the nucleus, it is not nuclear 
shadowing, since the conventional nuclear shadowing at low X2 does not require the invasion 
of an external proton to initiate the cascading interactions of gluons among the broken 
nucleons. We have no prediction on the nature of the effect of this type of depletion. If 
the experimental data on the production rate differ significantly from the curves in Fig. 12 
in the Xp < region, then there will be strong motivation to consider this unusual type of 
gluon depletion. 
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5 Nucleus-nucleus Collisions 



Having investigated the pA collision problem above, and finding the necessity to consider the 
depletion of gluons in the projectile before the production of cc states in the xf > region, 
it is natural to ask next what the implication would be for the J/ip suppression problem 
in AB collisions. Even without any detailed calculations it is straightforward to infer that 
there will be enhanced suppression at large xp. However, the currently available data on J/ip 
suppression in nuclear collisions are from CERN-SPS, and are limited to the central rapidity 
region. It is therefore our burden to show that the gluon depletion mechanism discussed in 
Sec. 3 is consistent with the existing nuclear collision data at small xp. 

In Sec. 3 we have used z mainly as the average penetration depth in the target nucleus 



where cc is produced, although in the evolution equation ([12]) z is the actual path length. 
Since in this section we need to average the production point over all impact parameters in 
the AB collisions, we now restore z to be actual path length of a gluon in a nucleus. Thus 
we allow it to vary from 1 to 12 fm. We use Eq. flTBD to determine g n (z) at all odd integer 
points in the range 1 < z < 12 (for the sake of simplicity in fitting), and then fit them at 
each such value of z by the formula (j32|) . The resulting values of a^z) are shown by the 
points in Fig. 13, they are in turn fitted by quadratic equations of the form 

a>i{z) = ko + b a z + b i2 z 2 . (36) 
The result is shown by the curves in Fig. 13. The corresponding coefficients are 
&io = 0.974, b n = -0.175, b 12 = 0.0068, 

b 20 = 0.0866, 6 2 i = 0.402, b 22 = -0.0202, (37) 
b 30 = -0.066, 6 31 = -0.226, b 32 = 0.0142. 



Using Eqs. (pq) and (p7|) in (|34|), we can evaluate g (xi, z) at all x\ and z. 

The cross section for the production of J ftp in AB collisions can be calculated in the 
standard way. We shall just write it down as follows (see, e.g., Refs [El. 0): 



crj/ip = J d 2 bd 2 sdz A dz B p A (s, z A )p B (s - b, z B ) 



• / — - — ~9a(x 1 , L b - z B ) g B (x 2 , L A - z A ) N(x t , A)N(x 2 , B) 
J X\ x 2 

. e --4PA(L A +z A )+p B (L B +z B )] & ^ c _( x ^ x ^ ^ ( 3g ) 

where we have included the z AyB dependences in the gluon distributions <?a,b(^i,2, L ByA — z b ,a) 
and the nuclear shadowing functions N(xi,A) and N(x 2 ,B). L AjB are the path length 
through A(B) at the distances s and s — b, respectively, from the centers of the nuclei, i.e. 



La — (Ra ~ s ) , L B = (H B - s-b Y' 2 . (39) 

The gluon distribution Qa^Xi, L b — z B ) is given by [see Eq. (|34]) 1 

3 

g A (xi,L B -z B ) =Y, a i( L B- z B ){l-xi) 4+i , (40) 
i=i 
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and similarly, for g B (x 2 , z a), the coefficients a,i(z) in Eq. fl34"| ) are replaced by La — za, which 
is the distance that a parton in B travels in A before the production of cc at za- The 
distance that a cc state travels in A is La + za, and za is integrated from — La to 
For the energies at CERN-SPS, \/s ^18 — 20 GeV, and for xp ~ 0, the hard cross section 
o"gg-*cc(^i) £2) restricts the gluon momentum fractions to x\ ~ £2 ~ Mj/^/y/s ~ 0.16. 
The survival probability is 

cAB _ .-.AB /_^(0) 
D J/^ — a J/4>/ a j/i> 

= N AB f d 2 bd 2 s I A dz A f * dz B W(b, s, z A , z B ) (41) 



where 



W(b, s, z A , z B ) = G A {x 1 ,L B - z B )G B {x 2 ,L A - z A ) N(xi,A)N(x 2 ,B) 

. e -°~a.{pA(LA+ZA)+PB(LB+ZB)} {4=2) 



and Nab is the same integral in Eq. (fHP but with s, Za, z b ) replaced by 1. Ga(%i, L b — 
z B ) is as given in Eq. (]35D except that 2; is replaced by L B — z B . Because of both the gluon 
enhancement at x\^ — 0.16 and the anti-shadowing, the absorption cross section a a now 
has to be somewhat larger than before 0. An overall agreement with all the pA and AB 
collision data, except the Pb-Pb case, can be achieved with the use of one value of a a = 8 
mb. The result is given in Table 2 for the various AB cases. Fig. 14 shows how those values 



Table 2: Survival probability for various AB collisions 



AB 


pp 


pC 


pAl 


pW 


pU 


OCu 


ou 


su 


PbPb 


S 


1 


0.85 


0.82 


0.70 


0.68 


0.67 


0.58 


0.56 


0.49 



compare with the experimental data |L4| by the straightline segments that connect those 
successive points. It is evident that apart from the Pb-Pb case the agreement with the data 
is satisfactory. Thus we can conclude that the gluon depletion mechanism used to treat the 
pA problem leads to no disagreement with the AB collisions — except that the Pb-Pb case 
remains as an anomaly. 

Finally, we compute the xf dependence of the J/ip suppression factor for just one nuclear- 
collision case as an example, which we take to be Pb-Pb. We consider two energies: = 
160 GeV and y/s = 60 GeV for RHIC. Except for the kinematics in <J 99 ^ cg that affects the 
values of X\ and x 2 , the cross section and survival factor can be calculated as before. The 
results are shown in Fig. 15, which exhibits a substantial degree of suppression at large 
xf- Any data at large xp would put considerable constraint on the models that attempt to 
explain the anomalous suppression at small xp- 
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6 Conclusion 



In our attempt to understand the enhanced suppression of J /if) at large xf, we have found 
that the depletion of gluons at large X\ in the projectile is the most natural explanation for 
the effect. We have proposed an evolution equation for the gluon distribution as the gluons 
propagate in a nuclear medium. The depletion at high x\ contributes to a mild growth of 
the gluon distribution at small x\. However, that growth does not lead to any contradiction 
with the existing data on J /if) suppression at mid-rapidity. Indeed, we have gone further to 
show where to find informative clues on the dynamics of suppression at large (positive and 
negative) Xp in both pA and AB collisions. 

What we have done here is only a modest first step towards understanding parton evo- 
lution in nuclear matter. While concentrating on the gluons, we have ignored the influence 
of the quark sector, a subject to be investigated at a later point. The depletion of quarks 
at large x reveals itself in the suppression of dilepton and leading meson production at large 
xf, the experimental evidences for which exist, though in subtle ways. Because of the con- 
servation of Fermion number, the degradation of the quark distribution at large x cannot 
be substantial. Nevertheless, the influence on the gluon distribution at small x may not be 
negligible. 

An important implication of this work is that in pA or AB collisions the concept of a nu- 
cleon propagating through nuclear matter as an identifiable, fixed entity needs modification. 
The usual notion that in nuclear collisions the total transverse energy E? is proportional to 
the number of nucleon-nucleon collisions would seem to have difficulty in reconciling with 
the insistence that the nucleons remain unaltered, if each inelastic collision of the nucleons 



contributes a fraction of their energies to Et- The wounded nucleon model |L5| makes a 
crude approximation of what goes downstream as an average quantity that is different from 
the incident nucleon, but ignores the way it changes as it propagates. Our evolution equa- 
tion indicates that the parton flux changes continuously and may emerge with a profile that 
cannot be identified with that of a free nucleon in any sensible comparison. The revelation 
made by this understanding will undoubtedly affect many aspects of high-energy nuclear 
collisions. 



Acknowledgment 

We are grateful to Kari Eskola for providing us with the Fortran codes for EKS98. This 
work was supported, in part, by the U.S. -Slovakia Science and Technology Program, the U. 
S. National Science Foundation under Grant No. INT-9319091 and by the U. S. Department 
of Energy under Grant No. DE-FG03-96ER40972. 



12 



References 

[1] R. C. Hwa, J. Pisut, N. Pisutova, Phys. Rev. Lett. 85, 4008 (2000). 

[2] R.C. Hwa, J. Pisut, and N. Pisutova, Phys. Rev. C 56, 432 (1997); 58, 434 (1998). 

[3] R. Vogt, Phys. Rep. 310, 197 (1999). 

[4] D. Kharzeev and H. Satz in Quark-Gluon Plasma 2, edited by R.C. Hwa (World Scien- 
tific, Singapore, 1995). 

[5] M. J. Leitch et al, (FNAL E866/NuSea Collaboration) Phys. Rev. Lett. 84, 3256 (2000). 
[6] Y.B. He, J. Hiifner, and B.Z. Kopeliovich, Phys. Lett. B 477, 93 (2000). 
[7] F. Arleo, P.-B. Gossiaux, T. Gousset and J. Aichelin, Phys. Rev. C 61, 054906 (2000). 
[8] R. Vogt, Phys. Rev. C 61, 035203 (2000). 

[9] M. J. Leitch et al. (FNAL E866/NuSea Collaboration) |nucl-ex/ 9909007 . 

[10] K.J. Eskola, V.J. Kolhinen, and PV. Ruuskanen, Nucl. Phys. B 535, 351 (1998). 

[11] K. J. Eskola, V. J. Kolhinen, and C. A. Salgado, Eur. Phys. J. C 9, 61 (1999). 

[12] Yu Dokshitzer, Sov. Phys. JETP 46, 1649 (1977); V. N. Gribov and L. N. Lipatov, Sov. 
Nucl. Phys. 15, 438, 675 (1972); G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298 
(1977). 

[13] R.C. Hwa, Phys. Rev. Lett. 52, 492 (1984). 

[14] M. Gonin et al, Nucl. Phys. A610, 404c (1996); M. C. Abreu et al. (NA50 Collabora- 
tion), Phys. Lett. B410, 337 (1997); B447, 28 (2000). 

[15] A. Bialas, M. Bleszynski, and W. Czyz, Nucl. Phys. Bill, 461 (1976). 



13 



Figure Captions 

Fig. 1 Nuclear shadowing and ant i- shadowing factor taken from EKS98 [|ll|] for Q 2 = 10 
GeV 2 . The curves are fits by the simple formula in Eq. (|9]). 

The function used to fit EKS98 data points at A = 100. 

K n calculated from the moments at discrete n and fitted by the formula in 
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curve 


Fig. 
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Fig. 
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Fig. 
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Fig. 
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The z dependences of K n (z). The lines are straight-line fits in the large z region. 
Q n as calculated from Eq. (|25|). 

The points of Q n are determined from Eq. (HB) or Fig. 5 at integer n values; the 



The kernal Q(y) as calculated from Eq. (|30] 



a(xp) for J/tj) production. The data points are from Ref. ||; the curves represent 
our result. 

Fig. 11 a(xp) for if)' production. The data points are from Ref. ||; the curves represent 
our result. 

Fig. 12 The survival probability for J/ip production in pW collisions at different energies 
for the entire range of Xp. 

Fig. 13 The z dependences of the coefficients a,i(z). 



Fig. 14 The survival probability in AB collisions. The data points are from Ref. [14 



the line is composed of straightline sections connecting the calculated points listed in 
Table 2. 

Fig. 15 The survival probability for J/ip in Pb-Pb collision for all Xp at two energies. 
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